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Abstract 

We apply a non-linear flow response formalism to the recently measured event plane correlations. 
We find that the observed event plane correlations can be understood as an average of the linear 
and quadratic response. 



1. Introduction 

The Quark-Gluon Plasma (QGP), as observed in heavy-ion collisions at RHIC and the LHC, 
exhibits strong collective flow JT] [2] [3] SJ . The observed flow pattern is characterized by long 
range correlations in the final state particle spectrum. Indeed, ATLAS recently measured event 
plane correlations |5|, and significant two-plane and three-plane correlations were observed. In 
comparison to the v„ measurements, the event plane correlations provide additional insight into 
the origin of initial state fluctuations and additional constraints on the shear viscosity of the 
QGP. Although these correlations can be simulated with event-by-event hydrodynamics [[6), it is 
important to explain the observed correlations without elaborate computer models. The aim of 
this paper is to explain the correlations, by taking into account the non-linear mixing of modes 
through quadratic order in "single-shot" hydrodynamics QUO. 



2. Methodology 

Event plane correlations can be obtained in single-shot hydrodynamics using a quadratic 
response formalism (9). As detailed in |9|, the response formalism uses the cumulant expansion 
(as opposed to a moment expansion) to characterize the initial geometry IfTUl . The medium 
response to a given cumulant is calculated at linear and quadratic order by deforming a Gaussian 
distribution and calculating the response coefficients with ideal and viscous hydrodynamics. 

For example, the fourth order cumulant C\ and angle <5>4 are closely related to the fourth 
order eccentricity ^r 4 e i4 ^y, and characterize the quadrangular deformation of energy density in 
the transverse plane 

C 4 e i4 ^ = (rV' 4 ^)-3(rV' % ) 2 . 

Using cumulants rather than moments to characterize the initial state geometry leads to different 
correlations between the angles O n . Fig. [1J shows the correlations between <t>2 and O4 with 
cumulant and moment based definitions. 
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Figure 1 : (a) Initial participant plane <I>2 and <I>4 correlation, as a function of centrality. Simulations are done with PHO- 
BOS Glauber model with <D 4 from both cumulants and moments definitions. The data points are extracted from ATLAS 
1 5 1 for comparison, (b) Event plane ¥2 and Y2 correlation, with predictions made by ideal and viscous hydrodynamics. 



The particle spectrum is expanded in harmonics 
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and each harmonic and angle, z„ = v n e~ m ^", will receive contributions from both the linear and 
quadratic response to a given cumulant. For v 4 for example, the dominant quadratic response 
is proportional to (e2e' 2 ° 2 ) 2 , and the linear response is proportional to the fourth order cumulant 

cV 4a>4 

Z4 = VA e-W* = w 4e - i4 ° 4 + w 4(22) e- i4 ^ , (2) 

where w 4 and w 4 (22) indicate linear response to C4 and the non-linear response to e^. 

It is straightforward then to express the event plane correlations using the complex form of 
harmonic flow z„ in Eq. ([TJ, e.g. 



cos 4(*P 4 - >P 2 ) = 



COS 4(04 — C*2)W4 + W4(22) 
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Eq. Q shows the typical structure of event plane correlations in the framework of non-linear 
response. In the numerator and denominator one term is determined by the linear response, w 4 , 
and one term is determined from quadratic response, w 4 (22)- In more involved cases, where higher 
order non-linear corrections are also included, there are interference terms between the linear and 
quadratic response. The combined effects of all these linear, non-linear, and interference terms 
determines the behavior of the event plane correlations. 

Fig. |TJb) shows that if linear response vv 4 dominates the v 4 signal, the event plane correla- 
tion reduces to its corresponding initial participant plane correlation, (cos4(0 4 - O2)}. On the 
other hand, if the v 4 signal is dominated by the non-linear response, w 4 (22), the event planes ¥4 
and *?2 would be perfectly correlated. Referring to these two limits as the linear and non-linear 
limits respectively, the observed event plane correlations can be understood as an average of 
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Figure 2: Two-plane correlations of event plane correlations measured at ATLAS 1 5 1. 
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Figure 3: Three-plane correlations of event plane correlations measured at ATLAS|5 |. The non-linear flow generations 
we considered in these calculations are: H>3(i2) for V3, w^zi) for V4, H"5(23) for V5, and ^6(33), W6(24. W6(222) for V6- Linear 
flow generation ivg is ignored because of its comparably small contributions. 
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these two limits, with the actual value controlled by the relative magnitudes of the two response 
coefficients, W4 and ^4(22). When these response coefficients are computed with ideal and vis- 
cous hydrodynamics and averaged as in Eq. (J3J, the event plane correlations seen in Fig. [TJb) 
are quantitatively reproduced by the response formalism. A more complete comparison of the 
response formalism to the two and three plane correlations is given in Fig.[2]and Fig. [3] 

The non-linear flow response dominates correlation for non-central collisions and for larger 
values of rj/s. Thus, we find that the observed non-trivial event plane correlations become 
stronger toward non-central collisions as seen in Fig.|2]and Fig. [3] Also as the shear viscosity is 
increased, the predicted correlations increase and approach the non-linear limit seen in Fig.[TJb). 
The observed correlation patterns involving <Df, have a rich structure. In the response formalism 
this rich pattern of correlations is due to the fac that two quadratric response coefficients, W6(24> 
and vt>6p3), contribute to Vg. 

3. Summary and discussion 

Using the PHOBOS MC-Glauber model, rj/s = l/4n, and a freeze-out temperature Tf Q = 
150 MeV, the non-linear response formalism qualitatively reproduces the observed two-plane 
and three-plane correlations, with one exception, ((cos(2 < f ( 2 - 6^3 + 4 v f , 4))). This correlator is 
currently under investigation. Since we have included only the quadratic response to the lowest 
deformations, this level of agreement is acceptable. In addition, a preliminary comparison of our 
predictions to event-by-event hydrodynamics also shows qualitative agreement [8 1. In particular, 
the response formalism reproduces the shear viscosity dependence of the correlations found in 
event-by-event hydrodynamics. A more complete comparison to event-by-event hydrodynamics 
and the data is reserved for future work. 
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